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ABSTRACT 

We show that a Magellanic Cloud origin for the warp of the Milky Way can explain most quantitative 
features of the outer HI layer recently identified by Levine, Blitz & Heiles (2005). We construct a 
model similar to that of Weinberg (1998) that produces distortions in the dark matter halo, and 
we calculate the combined effect of these dark-halo distortions and the direct tidal forcing by the 
Magellanic Clouds on the disk warp in the linear regime. The interaction of the dark matter halo 
with the disk and resonances between the orbit of the Clouds and the disk account for the large 
amplitudes observed for the vertical to = 0, 1,2 harmonics. The observations lead to six constraints 
on warp forcing mechanisms and our model reasonably approximates all six. The disk is shown to be 
very dynamic, constantly changing its shape as the Clouds proceed along their orbit. We discuss the 
challenges to MOND placed by the observations. 

Subject headings: Galaxy: structure, Galaxy: Disk, Galaxy: kinematics and dynamics 



1. INTRODUCTION 

The warp of the outer Milky Wa y, known since 1957 
ijKerr. Hindman fc Carpenter! Il957j) , has been quanti- 
tative ly determined for the first time by iLevine et al.l 
(2005). It can be described as a superposition of three 
and only three of the lowest order vertical harmonics 
of a disk: a dish-shaped to = 0, a integral- sign- shaped 
m = 1, and a saddle-shaped harmonic m = 2. The 
lines of nodes for each are close to coincident and nearly 
radial. The amplitudes of each reaches 7-10% of the 
radius of the disk. A number of possible warp produc- 
ing mechanisms have been suggested including long-lived 
eigenmodes, forcing by halo triaxiality, persistent cold- 
gas accretion, and tidal excitation. We show here that 
the origin of this warp can be well-described as the tidal 
interaction of the Magellanic Clouds with the disk and 
dark matter halo of the Milky Way. The interaction of 
the dark matter halo with the disk and resonances be- 
tween the orbit of the Clouds and the disk account for 
the large amplitudes of the three harmonics and their 
approximate shape and orientation. 

2. OBSERVATIONS 

ILevine et alJ 1)20051) found that a dynamical model for 
the warp must satisfy six observational constraints: 1) 
the three lowest-order harmonics: to = 0, 1 and 2 are 
necessary and sufficient to describe the global shape of 
the warp; higher-order global harmonics are typically an 
order of magnitude or more weaker; 2) the to = 1 warp 
has the largest amplitude everywhere in the outer disk; 3) 
the to = and to = 2 warps are comparable in amplitude 
to the m = 1 warp, but are smaller at all radii; 4) the to = 
1 warp has a measurable amplitude at the galactocentric 
radius of the Sun, R , but the to = and to = 2 warps 
begin near the edge of the stellar disk, at R — 2R ; 5) all 
three harmonics grow approximately linearly with radius, 
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reaching amplitudes of 1-2 kpc at about R — 30 kpc; and 
6) the amplitudes of each of the harmonics reach 5-10% 
of the radius of the disk. The lines of maximum descent 
of the to = 1 and to = 2 warps are coincident within 
about 12 degrees in galactocentric azimuth <f> and show 
little evidence of precession. The lines are located near 
<j) = 90°. 

3. METHODOLOGY 

We use the procedure described by iWeinberel l)1998l 
W98) to couple the halo response to the tida l excitation 
theory presented bv lHunter fc Toomrel (|1969t HT). This 
assumes that the disk remains thin and that gas dissi- 
pation is unimportant for the dynamics. HT derived an 
expression for the bending of a thin disk by the Magel- 
lanic Clouds and concluded that the direct excitation of 
the disk by the Clouds produces a warp of only a few 
hundred parsecs, an order of magnitude less than the 
amplitude of the obser ved warp. W98 w edded the halo 
excitation presented in Weinberg (1989) to the HT ap- 
proach, allowing the disk to feel both the tidal field from 
the Clouds directly as well as the force from the dark- 
matter halo wake excited by the Clouds. The assumption 
of linearity limits the predictions to modest amplitudes. 

We use perturbation theory rather than N-body simu- 
lation because of the intrinsic difficultly and subtlety in 
obtaining accurate multiple time scale results from a sim- 
ulation. The excitation hierarchy of satellite orbit— >halo 
wake— >disk bending modes results in multiple interleaved 
time scales: the orbital periods of the Clouds, the pat- 
tern speeds of the halo wakes, and the pattern speeds of 
the bending modes. We also have multiple spatial scales. 
For example, a simulated N-body disk must be capable 
of supporting the bending modes. In addition, particle 
simulation s can degrade the reson ant dynamics as de- 
scribed by iWeinberg fc Katd l|2005fl . With such difficul- 
ties, one should verify the dynamics of each component 
of the mechanism before putting everything together in 
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an N-body simulation, although this is rarely done. 
4. MODEL 

The perturbation theory includes the force from an 
extended satellite. However, the halo and disk excita- 
tion depend only on the lowest-order hamonics while the 
spatial extent changes the higher-order harmonics, and 
thus the extent plays little role. Therefore, as long as the 
Clouds remain bound, their masses may be added for our 
estimates. We use the disk pro file from W98. Our halo 
is an NFW profile with c = 15 l|Navarro et alJll997l) and 
virial mass 20 times the disk mass. The orbital plane was 
computed as described as in W 98. We used the radial ve- 
locity and proper motion frorn[Kallivavah^^yj2005) 
and the distance modulus from lFreedman et, alJ l)2001f ) to 
derive a space velocity and computed the resulting orbit 
in the sphe rical dark matter h alo. We adopted the LMC 
mass from iWesterhmdl l|1997t) of 2 x 1O 1O M , although 
there is considerable variance in the literature. The cur- 
rent position of the LMC is shown as is the current state 
of the warp. This orbit will carry it toward the NGP. 

Although our calculation assumes a collisionless and 
vertically thin medium, it is more generally applicable 
for the following reasons. Firstly, although the disper- 
sion relations for a multicomp onent and collisionless me - 
dia differ at small scales (e.g. Uod ll996: Rafikov 2001), 
the large-scale warp will be governed by inertia and the 
gravitational restoring force, not by local pressure. In 
addition, the disk self gravity is dominated by the inner 
stellar disk; the outer gas layer plays only a minor role 
in establishing the modes. Therefore, our collisionless re- 
sults are likely to be similar to a multicomponent calcula- 
tion. Secondly, the vertical restoring force depends very 
weakly on the thickening as long as the vertical degree 
of freedom does not couple to the bending (this may be 
demonstrated by straightforward but tedious algebra). 
Therefore, the modes will be largely unchanged by thick- 
ening over short time scales. However, the challenging 
problem of vertical coupling is important and remains to 
be investigated thoroughly. 

The warp is a very dynamic structure based on the 
temporal evolution of the model. This can be seen in 
the AVI file of the simulations which can be found at 
http : //www . astr o . umass . edu/^ weinberg/lmc. Also 
included at this site are comparisons of the m — 0,1,2 
evolution. Rather than a static structure that might be 
expected for a warp in response to a triaxial halo, a warp 
that results from the Magellanic Clouds is continuously 
changing shape because of the varying amplitudes and 
phases of the various modes. The image looks rather like 
a flag flapping in the breeze as the Clouds completes an 
orbit of the Milky Way. 

5. COMPARISON WITH OBSERVATIONS 

Figure ^ shows the shape of the warp for both the 
model and the data. We plot only m = 0, 1, 2 from the 
data, ignoring the weak but significant m = 10, 15 terms. 
The overall agreement is quite good, although there are 
some differences. The model does not, for example, have 
a minimum that is as extensive in Galactic azimuth as the 
data. Lines of maximum descent for m — 1, 2 in both the 
simulations and the HI data analysis are also shown in 
the figure. They are separated by about 20° in the model, 
a bit more than, but close to the 12° in iLevine et alJ 



( 2005) . Both sets of lines are oriented in the same sense: 
close to 4> — 90°. Neither set of lines shows evidence of 
significant variation with galactocentric radius. 

Figure shows the amplitudes of the first three har- 
monics, m = 0, 1, 2, in both the model and the data. We 
increased the Clouds' mass by 33% for a better fit. As in 
the data, the m — 1 in the simulations is the strongest, 
and increases nearly linearly out to the edge of the disk. 
At large R, where the inertia is sm all, the linear th eory is 
expected to overpredict the warp l)Tsuchival l2002) . The 
simulations also show a weak response of m = 0, 2 out 
to about 15 kpc, and then increasing nearly linearly, but 
with approximately the same amplitude for both, pro- 
viding a reasonable representation of the data. 

The amplitude ordering of these vertical harmonics 
has a natural physical explanation. The force from 
the halo wake and the satellite may described by three- 
dimensional harmonics and these affect the warp height 
as follows. A spherically symmetric m = halo distor- 
tion does nothing to the warp. An m = 1 halo distortion 
will tend to accelerate the disk in the vertical direction; 
the differential acceleration of the disk results in a ver- 
tical m — "dishing" . An m = 2 halo distortion will 
attract the disk upward and downward in a reflection- 
symmetric way causing the classic integral-sign warp. 
Higher-order symmetries may be deduced by similar ge- 
ometric considerations. The power in the halo excitation 
drops off as an inverse power of the harmonic order and 
only the lowest order terms have features well inside the 
satellite orbit. Conversely, the existence of these higher- 
order harmonics with a power-law drop off is a natural 
consequence of this tidal theory and is consistent with 
the data. 

We fixed the disk and the halo mass inside of the virial 
radius while adjusting the satellite orbit and halo concen- 
tration and found the following trends. First, the halo 
wake and its pattern speed are determined by the halo 
concentration. The disk bending modes have a natural 
set of frequencies for a given halo. These will be max- 
imally excited when forced by the halo wake at or near 
harmonics of this natural frequency. The ratio of m = 2 
amplitude to the m = 1 amplitude is maximized for a 
halo concentration c « 10. For an NFW profile, this 
puts r s w 30 kpc. This is very close to the ACDM es- 
timates for the Milky Way concentration. Secondly, the 
orientation of the response depends on whether the near- 
est resonance is larger or smaller than natural frequency. 
Therefore changing the satellite orbit, which changes the 
forcing frequency, affects both the amplitude of the am- 
plitude and orientation of the warp response. In short, 
the warp depends on a "clockwork" of frequency rela- 
tionships which depends on the satellite orbit, the dark- 
matter halo and the disk. A pericenter larger than the 
current 49 kpc estimate shifts the position angle (PA) so 
that the warp peaks closer to PA=180 degrees. Similarly, 
smaller pericenter increases the amplitude also shifts the 
PA. We conclude that our model "prefers" our current 
fiducial MC model. We are not claiming that our fiducial 
model is the most probable amongst the distributions of 
allowed values but that a plausible choice of parameters 
corresponds to many features of the observed data. 
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Fig. 2. — The amplitudes of the m = 0,1,2 harmonics as a 
function of Galactic radius fro m the simulations (c urves) and from 
the HI Milky Way data (points ITevine et all2003) . 

6. DISCUSSION 
6.1. Comparison with N-body Simulations 

The predictions of W98 wer e checked by several groups 
using N-body simulations. IGarcia-Ruiz et all (2002) 
used a hybrid N-body particle-ring code and did not 
find the warp predicted in W98. They offer the failure of 
the linear theory as the culpr it but did not invest i gate a 
variety of models. Similarly, IMastropietro et all ( 2005) 
performed a simulation including both the gaseous, stel- 
lar and dark components of the Large Cloud and the 
Milky Way and remark th at the effe c t on t he Milky Way 
is negligible. Conversely, ITsuchival l)2002|) . using a hy- 
brid code that includes both a potential expansion and 
a tree code, obtained amplitudes of m = 0,1,2 warps 
in the larger halo model that show very good agreement 
with the observed amplitudes in the Milky Way. He also 
shows that the amplitude of the warp depends strongly 
on the dark halo mass model. 

It is difficult to reconcile these contradictory findings 
but three possibilities obviously occur: 1) the linear the- 



ory does not apply; 2) the simulations do not apply 
due to numerical difficulties; and 3) the two simulations 
with null warps have chosen unlu cky sets of parame- 
ters. Both lGarcia-Ruiz et alJ l|2002fl and lTsuchival (2002) 
choose methods that exp licitly treat the multiple scales. 
IGarcia-Ruiz et all l)2002|) used tilted rings to represent 
the disk and a tree code to represent the halo. Although 
the use of rings limits the investigation to the m = 1 re- 
sponse, it should be sufficiently sensitive to th e halo ex- 
citatio n without strong particle number issues. ITsuchival 
( 2002) used an expansion algorithm to represent the halo 
gravitational field and a tree code to repr esent the disk to 
better represent the multiple scales (see ITsuchival 12001 
for additional discussion of the dynamic range problems 
that motivate this choice). We feel that the discrepancy 
between the results of these two groups is most like Item 
(3): the Garcia- Ruiz et al. model is not particularly 
warp producing. The good quali tative correspondence 
between W98 and lTsuchival l|2002j) suggests that the lin- 
ear theory captures the underlying physics although is 
likely to differ in detail. For example, the linear theory 
over-predicts the height in the edge of the disk (Fig. EJ, 
and this motiv ates our truncation of the predictions in 
the outer disk. IMastropietro et all l)2005fl primarily em- 
phasize the affect of the Milky Way tides on the Large 
Cloud and do not tailor their approach to treat the mul- 
tiple scale problem per se. We suggest that their report 
of no disk warp results from Items (2) and (3). 

6.2. Other Explanations 

Most warp theories depend on the bending response of 
the disk and this response is strongly affected by the 
existence of a self-gravitating dark matter halo. The 
theories may be roughly grouped as follows: 1) bending 
modes may be primordial ly excited and persistent (e.g. 
Soar ke fc Casertanoll988(k 2) the response of the disk to 
the non-axisymmetric shape of the halo; 3) tidal excita- 
tion (as we have discussed here) and 4) the warp may be 
produced by the response of the disk to cosmic infall (e.g. 
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Jiang fe Binneviri999() . The t opic was nicely reviewed by 
Binnevll|199lir Siibseaiientlv. lNelson fc Tremain'el l|1995ft 
argued against long-lived modes. The triaxialit y of the 
Milky Way halo remains unc ertain. iHelmil (l2004f) re ports 
a prolate halo (q=1.25) while l.Tohnston et alJ (2005) finds 
an oblate halo (q=0.8-0.9) using Sgr dwarf constraints. 
Such modest triaxiality seems unlikely to produce the ob- 
served vertical m — 2 feature. In addition, cosmic infall 
more naturally produces tilted rings, an m = 1 feature; 
the observed vertical m = 2 may require a conspiracy of 
several inflow directions. 

6.3. Relevance for Modified Gravity 

Although the success of our model favors the existence 
of a dark-matter halo in Nature, many find it seduc- 
tive to modify gravity to produce the observed rota- 
tion velocities in the Galaxy without dark matter (mod- 
ified Newtonian dynamics o r MOND iMilgrorr] 1 1 983albt 
iBekenstein fc Milgromll984|) . Might the tidal theory also 
apply in MOND? It is beyond the scope of this letter to 
repeat our calculations using MOND, it seems plausible 
that direct forcing of the Galactic disk by the Clouds in 
MOND may provide warp amplitudes in excess of the 
original predictions without a dark halo (HT). Similarly, 
we would expect the disk modes and frequencies to be 
qualitively similar, the excess restoring force of the halo 
being produced by the MOND force. However, MOND 
would have to (1) admit bending modes with similar mor- 
phology to those in the Newtonian theory and (2) these 
modes would have to conspire to frequencies that couple 
them to the direct forcing by the Clouds in such a way 
that they assumed the same orientation as in the "clock- 
work" described in §5. Because this clockwork and the 
m = 2 to m = 1 amplitude ratio depends on the si- 
multaneous halo and disk excitation, agreement with the 
observations described in §2 seems rather unlikely and 
thereby disfavors MOND. We encourage detailed predic- 
tions. 

6.4. Relevance for other Galactic Systems 

A large fraction of other warped galaxies also show 
warps in their HI layers and a significant fraction of these 
are asymmetric. Warps in these systems are generally 
analy zed with a program such at ROTCUR (Begeman 
1989) or one of its derivatives, but invariably are forced 
to fit the to = 1 wa rp only. This paper, and the work of 
iLevine et alJ l)2005l) show that at least three harmonics 
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ought to be fit to the warps of external galaxies, espe- 
cially those with asymmetric warps, which can be caused 
by a superposition of these harmonics. The degree to 
which various harmonics are present in a warp can pro- 
duce important constraints on whether the warp is due 
to a satellite, a triaxial halo, cold gas inflow, or some 
primordial excitation. All warps may not be alike. 

6.5. Summary 

We have demonstrated a plausible mechanism for the 
excitation of the warp that explains all of its general fea- 
tures. A number of warp producing mechanisms have 
been explored elsewhere in addition to tides. It is possi- 
ble that different mechanisms may act in different galax- 
ies and possibly in concert in a single galaxy. Nonethe- 
less, our simple model nicely reproduces the observed 
features in the Milky Way HI gas layer. The existence of 
massive companions, the Magellanic Clouds, and the pre- 
diction from linear perturbation theory and at least one 
corroborating N-body simulation suggests that a tidal 
explanation is viable. Our model depends on the gravi- 
tational response of the halo and thereby suggests that 
the dark matter is not an artifact of modifying the laws 
of gravity. Conversely, given a dark halo, we argue that 
the tide from the Magellanic Clouds must be affecting 
the Milky Way disk, and, given the quality of the agree- 
ment, it seems to be the dominant mechanism. This 
model then promises an additional constraint on the dis- 
tribution of dark matter. Although we have emphasized 
the gas layer response beyond the stellar disk, the effect 
of other satellite encounters such as the recent accretion 
of Sagittarius dwarf may be detectable in future high- 
resolution surveys and may help determine the proper- 
ties of the inner halo. Warp observations are important 
because promise to reveal aspects of the dark matter dis- 
tribution that are otherwise observationally inaccessible. 
Further analysis of warps may provide more precise con- 
straints on the profile of dark matter in the Milky Way 
and nearby external galaxies. 
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